Is a co-rotating Dark Disk a threat to Dark Matter Directional Detection ? 
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Recent N-Body simulations are in favor of the presence of a co-rotating Dark Disk that might 
contribute significantly (10%-50%) to the local Dark Matter density. Such substructure could have 
dramatic effect on directional detection. Indeed, in the case of a null lag velocity, one expects an 
isotropic WIMP velocity distribution arising from the Dark Disk contribution, which might weaken 
the strong angular signature expected in directional detection. For a wide range of Dark Disk 
parameters, we evaluate in this Letter the effect of such dark component on the discovery potential 
of upcoming directional detectors. As a conclusion of our study, using only the angular distribution 
of nuclear recoils, we show that Dark Disk models as suggested by recent N-Body simulations will not 
affect significantly the Dark Matter reach of directional detection, even in extreme configurations. 

PACS numbers: 95.35.+d, 14.80.-j 



I. INTRODUCTION 



Within the standard Dark Matter halo paradigm, 
the local Dark Matter distribution is assumed to be 
smoothly spatially distributed and to be well-described 
by a Maxwellian velocity distribution. However, the hi- 
erarchical structure formation model indicates that the 
Galactic Dark Matter halo results from sucessive small 
halo accretions, thus directly linking its structure to its 
merging history. The presence of substructures in the 
Milky Way halo is inferred from recent results of N-body 
simulation PHI]. Such substructures may be classified 
as follows : Dark Matter tidal streams (spatially local- 
ized), debris flows (spatially homogenized but with ve- 
locity substructures) and a Dark Disk. The latter has 
received much interest since late sub-halo merging is ex- 
pected to lead to the formation of a co-rotating Dark Disk 
[g-Q that may affect the expected WIMP signal both in 
direct and directional detection. While the influence of 
the dark disk on Dark Matter signals has been exhaus- 
tively investigated for direct @, Il0r{l2j and indirect [l3[ 
detection, it is still unclear how it may affect directional 
detection. Following a previous work from A. M. Green 
fl2j , we aim at evaluating the influence of the presence 
of a co-rotating Dark Disk on the discovery potential of 
a forthcoming directional detector. In particular, we are 
interested in determining the values of the Dark Disk 
parameters for it to significantly affect the Dark Mat- 
ter reach of directional detection. In order to be model 
independent from the background energy modelling, the 
study has been done by considering only the angular dis- 
tribution of nuclear recoils dR/dU r . 
Since the pionner paper of D. N. Spergel [1J|, the contri- 
bution of directional detection to the field of Dark Matter 
has been adressed through a wealth of studies [1, [l5l - l36| . 
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Depending on the unknown WIMP-nucleon cross section, 
directional detection may be used to : exclude Dark Mat- 
ter [H, [r| , reject the isotropy hypothesis fl7l - [23| , dis- 
cover galactic Dark Matter with a high significance 0- 
I2H ] or constrain WIMP and halo properties [27l - r29j . For 
neutralino Dark Matter, a large fraction MSSM configu- 
rations with a neutralino lighter than 200 GeV/c 2 would 
lead to a significance greater than 3a (90% CL) in a 30 
kg. year CF4 directional detector 30]. 
In the following, we focus on the effect of a co-rotating 
Dark Disk on the potential of forthcoming directional 
detectors to discover Dark Matter 26]. The paper is 
organized as follows. Section [TT] presents the current 
knowledge on the Dark Disk. In particular, we define 
its parametrization used throughout. The directional 
framework is recalled in sec. IIII1 with emphasize on the 
directional statistic methods used to exploit forthcoming 
data. Then, for a wide range of Dark Disk parameters, 
we evaluate in IIVI the effect of such substructure on the 
discovery potential of upcoming directional detectors. 



II. A DARK DISK IN THE MILKY WAY 

Recent results from N-body simulation of Milky Way 
type galaxies have shown that merging satellite galax- 
ies may get dragged into the plane of their host galaxy 
This leads to a Dark Matter overdensity roughly 
matching the baryonic disk of the host galaxy and usu- 
ally in co- rotation with the latter 0, [1| . This Dark Mat- 
ter component is usually referred to as Dark Disk (DD). 
Hitherto, there is no observational evidence in favor of a 
Dark Disk in the Milky Way. In fact, a Dark Disk den- 
sity greater than 50% of the halo one (pdd/ph > 0.5) 
is excluded by an estimation of the dynamical surface 
mass density in the solar neighbourhood obtained from 
the kinematics of a sample of 1200 red giant stars [371 ]. 
The Dark Disk is generally considered as a cold substruc- 
ture for which the velocity distribution is decribed by an 
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isotropic Maxwellian distribution @- In such context, 
the astrophysical parameters relevant to the description 
of the Dark Disk are the density pdd, its co- rotational 
velocity Vdd and its velocity dispersion given by odd- 
The range of interest of these parameters must be inferred 
from the results of N-body simulations and compared to 
astrophysical constraints. For instance, F. S. Ling et al. 
@ have extracted a Milky way type galaxy from the re- 
sults of the RAMSES simulation Q. The velocity dis- 
tribution of Dark Matter particles within a 7 < R < 9 
kpc and \Z\ < 1 kpc is shown to be well fitted by a 
double gaussian along the <fi direction. The first one, 
corresponding to the Dark Matter halo component, is 
described by a null average speed and a velocity disper- 
sion Ohaio — 180 km/s. The second component, the Dark 
Disk, is described by a co-rotation velocity Vdd — 150 
km/s and a velocity dispersion odd — 85 km/s. Note 
that this description is in good agreement with 043- 
However, as there is no clear observational constraints 
on these parameters, we allow for a wide range in order 
to investigate the effect of a Dark Disk component on di- 
rectional detection. Unless otherwise stated, we consider 
hereafter the following Dark Disk parameter ranges: 



< PDD/ph < 1 
km.s -1 < Vdd < 220 km.s" 1 
7 km.s -1 < (Tdd < 155 km.s" 1 



III. 



DIRECTIONAL DETECTION 
FRAMEWORK 



(1) 



There is a worldwide effort toward the development 
of a large TPC (Time Projection Chamber) devoted to 
directional detection j38l - l42| . In the following, we ex- 
emplify our Dark Disk study considering a MIMAC-like 
detector corresponding to a low exposure (30 kg. year) 
CF4 TPC allowing three dimensional track measure- 
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The two dimensional directional recoil rate 
d 2 R/dE r dVL r is given by [35[ : 



d 2 R 

dE r dfl r 



Ai:m x m' 2 . 



F 2 (E r )f(v min ,q), 



(2) 



with m x the WIMP mass, m r the WIMP-nucleus 
reduced mass, po the local Dark Matter density, <7q 
the WIMP-nucleus elastic scattering cross section, 
F{E r ) the form factor (using the axial expression 
from [43]), u m in the minimal WIMP velocity required 
to produce a nuclear recoil of energy E r and q the 
direction of the recoil momentum. Finally, f(v m i ni q) is 
the three-dimensional Radon transform of the WIMP 
velocity distribution f(v). As the Radon transform 
is a linear application, one can simply add the host 
halo and the Dark Disk contribution to the directional 
event rate. The angular distribution dR/dfl r is thus 



obtained by integrating the double-differential spec- 
trum over the energy range chosen to be E r = [5, 50] keV. 

There are several approaches to exploit the forthcom- 
ing directional data. Either the data analysis may aim 
at rejecting the isotropy hypothesis [171423T ] or at discov- 
ering galactic Dark Matter with a high significance (2(| 
via a profile likelihood ratio test. In the following, we 
investigate the effect of a Dark Disk contribution on the 
expected significance of a discovery of Dark Matter with 
directional detection. To do so, we use a twofold ap- 
proach, following [20L [26j. 

First, we use a profile likelihood ratio test statistic, as 
presented in [26| and brielfy recalled hereafter for the 
reader's convenience. In order to remain independent 
from the background energy spectrum modelling, only 
the directional information is considered in the follow- 
ing, i.e. the angular distribution of the recoiling events 
dR/dil r . Noting a p the WIMP-proton cross section and 
Rb the background rate, the likelihood function is given 
by, 

J?(<r p ,R b )= {fXs+ N f N e-^+^ 
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where pb — Rb X £ and p s corresponds to the number 
of expected background and WIMP events respectively, 
where £ corresponds to the exposure. iV is the total num- 
ber of observed events, R n refers to the direction of each 
event while the functions S and B are the directional 
event rate dR/dfl r of the WIMP and background events 
respectively. Following recent studies on the angular dis- 
tribution of muon induced neutrons (44| , the background 
angular distribution B is assumed to be isotropic in the 
galactic rest frame. Note that, contrary to a previous 
work [2^ | the astrophysical uncertainties are not taken 
into account in the estimation of the significance to allow 
fair comparison between the two statistical approaches. 
Only the background rate is taken as a nuisance param- 
eter. 

In a frequentist approach, the significance of a new pro- 
cess is commonly estimated by using the profile likelihood 
ratio test [45[ . It corresponds to a hypothesis test of the 
null hypothesis Hq (background only) against the alter- 
native Hi which includes both background and signal. 
As discussed in [45[ the test statistic in the case of a 
discovery is defined as follows: 
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V DD [km/s] V DD [km/s] 

FIG. 1: Mean significance E(Z) as a function of Vdd, the rotation velocity of the Dark Disk at Solar radius. Left : the result 
is presented for various values relative density poo/pH and a fixed velocity dispersion ctdd = 85 km/s. Right : the result 
is presented for various values of the velocity dispersion a DD and a fixed value of the relative density pdd/ph = 0.5. These 
studies has been done for 50 GeV/c 2 WIMP mass and a fixed value of 100 WIMP events (from the halo and the Dark Disk). 



Hence, a large value of qg implies a large discrepancy 
between the two hypothesis which is in favor of a dis- 
covery (Hi). As f(qo | Hq) follows a \\ distribution, the 
discovery significance Z is simply defined as Z = \J q ohs , 
in units of a [45j |. 

The second approach, first introduced by B. Morgan 
et al. [20] , is based on a generic test of isotropy following 
the mean recoil deviation (cos 9) such as: 

1 N 

(cos 9) = — cos 0» (6) 

i=l 

where 6{ is the i th angle between the recoil and the 
Cygnus direction, and N is the number of measured re- 
coils. Note that this test is, by definition, coordinate 
system dependent, such as the previous one as the main 
recoil direction (£,b) (see (25[) is not considered here as 
a fitting parameter. 

Eventually, one can evaluate the significance of an ob- 
served anisotropy by computing the distributions of 
(cos 9) for both Hq corresponding to the background 
(isotropic) and Hi the alternative. It is worth noticing 
that the use of the variable (cos 9) is particularly interest- 
ing in the case of directional detection of Dark Matter as 
the expected signal should exhibit a dipole feature hence 
maximizing the deviation between Hq and Hi. 

IV. INFLUENCE OF A CO-ROTATING DARK 
DISK 

In order to investigate the effect of a Dark Disk com- 
ponent on the expected significance of a directional dark 
matter detection, we allow for a wide range on the Dark 
Disk parameters, see eq. [TJ and we evaluate, for each 



configuration, the expected significance for a 30 kg. year 
MIMAC-like detector. We highlight the fact that for a 
co-rotating Dark Disk to contribute to the data, the en- 
ergy threshold must be low and/or the WIMP mass large. 
For concreteness, we present a case study for a 50 GeV/c 2 
WIMP mass and a total of 100 WIMP events. Figure Q] 
(left) presents the mean significance E(Z) as a function 
of Vd d , the rotation velocity of the Dark Disk at Solar 
radius. The black dashed line corresponds to the no Dark 
Disk case. The result is then presented for various values 
of the relative density pdd/ph- The general feature is 
that the mean significance is decreasing when increasing 
the co-rotating velocity of the Dark Disk at Solar radius 
as it results in a loss of directionality. This effect is even 
stronger when increasing the Dark Disk contribution, i. e. 
for large values of the relative velocity pdd/ph- Inter- 
estingly, a co-rotating Dark Disk can boost the mean sig- 
nificance of a Directional Dark Matter detection. Indeed, 
for a velocity dispersion udd = 85 km/s and a WIMP 
mass of 50 GeV/c 2 , one can see that for rotation velocity 
Vdd < 140 km/s and for any relative density, the mean 
significance obtained is greater than the one obtained in 
the no Dark Disk case. This enhancement of the signifi- 
cance at low rotation velocities can be explained by the 
fact that the Dark Disk is a colder structure than the host 
halo, i.e. has a smaller velocity dispersion, implying an 
even more anisotropic recoil angular distribution. Hence, 
a co-rotating Dark Disk will not necessarily degrade the 
expected performance of directional detection. Of course, 
for a perfectly co- rotating Dark Disk (Vdd — vq = 220 
km/s), whatever the velocity dispersion, the recoil angu- 
lar distribution induced by the Dark Disk is necessarily 
isotropic. Only the contribution of the Dark Disk to the 
total number of events will change. 

Figure [T] (right) presents the mean significance as a func- 
tion of Vdd- F° r an y value of the velocity dispersion 
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FIG. 2: Relative asymmetry s# in the plane (Vdd,(t v ) for a relative density pdd/ph = 1/3 (left) and pdd/ph = 1 (right). 
The solid and dashed lines correspond to isocontours of a relative asymmetry equal to ±1% and ±5% respectively. These 
studies have been done for 50 GeV/c 2 WIMP mass. 



<JDDi the mean significance is continuously decreasing 
with the rotation velocity of the Dark Disk. However, 
the case odd — 35 km/s (red solid line) tends to the no 
Dark Disk limit due to the fact that the contribution to 
the total number of WIMP events from the Dark Disk 
falls quickly to zero for Vdd > 140 km/s. This also ex- 
plains the rapid decrease of the signficance enhancement 
in the range — 100 km/s. For larger velocity dispersions, 
the mean significance does not tend to the no Dark Disk 
limit as the contribution of the Dark Disk to the total 
number of events remains non negligible. Interestingly, 
one may note that the range of the values of Vd d induc- 
ing an enhancement of the significance depends strongly 
on odd- Indeed, for large values of the velocity disper- 
sion, the Dark Matter signal gets closer to an isotropic 
distribution. This constatation implies that lower is the 
velocity dispersion, larger is the range in values of Vdd 
allowing for a boost of the directional signature. As a 
conclusion, larger is the velocity dispersion of the Dark 
Disk, weaker is the directional discovery significance, ex- 
cept for the case of odd = 35 km/s as discussed above. 
However, note that the value of <tdd = 141 km/s is ex- 
tremely large with respect to the recent results from N- 
Body simulations. Hence, for a 50 GeV/c 2 WIMP mass, 
one could expect that a co-rotating Dark Disk could have 
a positive, though small (~ 10%), effect on the directional 
detection of Dark Matter. 

For completeness, we studied the evolution of the mod- 
ifications of the angular distribution dR/d£L r for various 
Dark Disk parameter values. For this purpose, we defined 
the relative asymmetry si as: 

(cos6>) - (cos 6} H 
* = ^Th (7) 

where {cos9)h corresponds to the mean recoil deviation 
obtained in the no Dark Disk case (pdd — 0) . Note that 
considering the standard halo model and a WIMP of 50 



GeV/c 2 , we found {cos0)h — 0.51. Figure [21 presents 
the relative asymmetry si in the plane (Vdd,&® D ) for 
a relative density pdd/ Ph = 1/3 (left) and pdd/ph = 1 
(right). One may notice that there are three different re- 
gions : no effect (the 1% region), a directional discovery 
enhancement region (low Vdd) and a region for which 
the Dark Disk weakens the directional signature (high 
Vdd together with a high odd value). The relative 
density only affects the amplitude of si ', note that the 
latter spans the range [-18,10] for pdd/ph = 1/3 and 
[-40,23] for pdd I Ph = 1- This also affects the area of the 
no effect region which decreases with increasing value of 
Pdd I ' Ph- Interestingly, one can notice that most of the 
Dark Disk models suggested by N-Body simulations lie 
in the no effect region. Only extreme, yet unrealistic, 
Dark Disk models may affect significantly the directional 
signature. 

We evaluate the effect of the Dark Disk contribution to 
the Dark Matter reach of upcoming directional detectors. 
Following 26], we compute the directional reach in the 
(to x , cr p ) plane, i.e the lower bound of the 3tr discovery re- 
gion at 90% CL for the two approaches: profile likelihood 
(red lines) and mean recoil deviation (blue lines) . Figure 
[3]presents the discovery limit in the {m x , log 10 (<7p)) plane 
corresponding to two Dark Matter models: standard halo 
model only (solid lines) and and with an extreme Dark 
Disk model contribution {pdd/ Ph = 1, Vdd — 220 km/s, 
cdd = 106 km/s} (dashed lines). 

The conclusion of this study is twofold. First, we found 
that for both statistical tests, the effect of an extreme 
Dark Disk is only mild. Indeed, the directional reach is 
only degraded by a factor of 3 at high WIMP masses 
and not affected for light WIMP. Second, we found that 
the two statistical tests give similar results with a max- 
imal deviation of a few percent. However, it is worth 
emphasizing that there intepretation differ as the profile 
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FIG. 3: Discovery limit in the (m x , log 10 (<r p )) plane corre- 
sponding to two Dark Matter models: standard halo model 
only (solid lines) and and with an extreme Dark Disk model 
contribution {pdd/ph = 1, Vdd = 220 km/s, odd = 106 
km/s} (dashed lines). We compute the directional reach, i.e 
the lower bound of the 3a discovery region at 90% CL, for the 
two approaches: profile likelihood (red lines) and mean recoil 
deviation (blue lines). 



likelihood method favors the background plus signal hy- 
pothesis (H\ ) whereas the mean recoil deviation method 
rejects the isotropy hypothesis. 



V. CONCLUSION 

A co-rotating Dark Disk, as predicted by recent N- 
Body simulations, might contribute (10%-50%) to the 
local Dark Matter density, with a potentially dramatic 
effect on directional detection. In this letter, we have 
evaluated the effect of Dark Disk model on the discovery 
potential of upcoming directional detectors. We conclude 
that, if a co-rotating Dark Disk is present in our Galaxy 
and has the properties predicted by N-Body simulations 
Q , the discovery potential of directional detection would 
be strictly unchanged. Only an extreme, yet unrealis- 
tic, Dark Disk model might affect significantly the Dark 
Matter reach of upcoming directional detectors, by in- 
creasing the discovery limit by a factor of three at high 
WIMP mass (m x ~ 1000 GeV/c 2 ). 

Note that even if the impact of Dark Disk contribution 
to the local Dark Matter distribution only mildly affects 
the discovery potential of directional detection, it may 
significantly affect the mass and cross section determi- 
nation [27j . Indeed, as explained in [l2j |. WIMP events 
arising from the Dark Disk contribution will induce an 
excess at low recoil energies which can lower the estima- 
tion of the WIMP mass when considering a standard halo 
model. This highlights the fact that even if a co-rotating 
Dark Disk is not a threat to the discovery potential of 
directional detection, it has to be characterized in order 
to consistently constrain the Dark Matter properties. 



R. Teyssier, Astron. Astrophys. 385 (2002) 337 [18 

M. Lisanti and D. N. Spergel. larXiv:1105.4166l 

M. Vogelsberger et al, Mon. Not. Roy. Astron. Soc. 395 [19 

(2009) 797 

M. Kuhlen, M. Lisanti and D. N. Spergel, [20 
larXiv:1202.0007l [astro-ph.GA] . 

J. Diemand et al., Nature 454 (2008) 735 [21 

F. S. Ling, E. Nezri, E. Athanassoula and R. Teyssier, 
JCAP 1002 (2010) 012 [22 
J. I. Read, G. Lake, O. Agertz and V. P. Debattista, 
MNRAS 389 (2008) 1041-1057 [23 
J. I. Read, L. Mayer, A. M. Brooks, F. Governato and 

G. Lake, MNRAS 397 (2009) 44 [24 

C. W. Purcell, J. S. Bullock and M. Kaplinghat, Astro- 
phys. J. 703 (2009) 2275 [25 
T. Bruch, A. H. G. Peter, J. Read, L. Baudis and G. Lake, 
Phys. Lett. B 674 (2009) 250 [26 
F. S. Ling, Phys. Rev. D 82 (2010) 023534 
A. M. Green, JCAP 1010 (2010) 034 [27 
T. Bruch, J. Read, L. Baudis and G. Lake, Astrophys. J. 

696 (2009) 920 [28 

D. N. Spergel, Phys. Rev. D 37 (1988) 1353 [29 
J. Billard, F. Mayet and D. Santos, Phys. Rev. D 82 

(2010) 055011 [30 
S. Henderson, J. Monroe and P. Fisher, Phys. Rev. D 78 
(2008) 015020 [31 
C. J. Copi and L. M. Krauss,, Phys. Rev. D 63 (2001) 
043507 [32 



A. M. Green and B. Morgan, Phys. Rev. D 77 (2008) 
027303 

C. J. Copi and L. M. Krauss, Phys. Lett. B 461 (1999) 
43 

B. Morgan, A. M. Green and N. J. C. Spooner, Phys. 
Rev. D 71 (2005) 103507 

B. Morgan and A. M. Green, Phys. Rev. D 72 (2005) 
123501 

A. M. Green and B. Morgan, Astroparticle Physics 27 
(2007) 142 

C. J. Copi, L. M. Krauss, D. Simmons-Duffin and 
S. R. Stroiney, Phys. Rev. D 75 (2007) 023514, 

A. M. Green and B. Morgan, Phys. Rev. D 81 (2010) 
061301 

J. Billard, F. Mayet, J. F. Macias-Perez and D. Santos, 
Phys. Lett. B 691 (2010) 156-162 

J. Billard, F. Mayet and D. Santos, Phys. Rev. D 85 
(2012) 035006 

J. Billard, F. Mayet and D. Santos, Phys. Rev. D 83 
(2011) 075002 

S. K. Lee and A. H. G. Peter, JCAP 1204 (2012) 029 

D. S. M. Alve s, S. E. Hedri and J. G. Wacker, 
larXiv: 1204.54871 [astro-ph.GA] . 

D. Albornoz Vasquez, G. Belanger, J. Billard and F. 
Mayet, Phys. Rev. D 85 (2012) 055023 
N. Bozorgnia, G. B. Gelmini and P. Gondolo, 
larXiv:1111.636T1 [astro-ph.CO]. 

M. S. Alenazi and P. Gondolo, Phys. Rev. D 77 (2008) 



6 



043532 

[33] R. J. Creswick, S. Nussinov and F. T. Avignone, As- 

tropart. Phys. 35 (2011) 62 
[34] N. Bozorgnia, G. B. Gelmini and P. Gondolo, 

larXiv:1205.2333l [astro-ph.CO] . 
[35] P. Gondolo, Phys. Rev. D 66 (2002) 103513 
[36] M. Lisanti and J. G. Wacker, Phys. Rev. D 81 (2010) 

096005 

[37] C. M. Bidin, G. Carraro, R. A. Mendez and W. F. van 

Altena, Astrophys. J. 724 (2010) L122 
[38] S. Ahlen et ai, Int. J. Mod. Phys. A 25 (2010) 1 



[39] K. Miuchi et ai, Phys. Lett. B 686 (2010) 11 
[40] S. Ahlen et ai, Phys. Lett. B695 (2011) 124-129 
[41] D. Santos et ai, EAS Publications Series 53 (2012) 25-31 
[42] E. Daw et ai, Astropart.Phys. 35 (2012) 397-401 
[43] J. D. Lewin and P. F. Smith, Astropart. Phys. 6 (1996) 
87. 

[44] D. Mei and A. Hime, Phys. Rev. D 73 (2006) 053004 
[45] G. Cowan, K. Cranmer, E. Gross, O. Vitells, Eur. Phys. 
J. C71 (2011) 1554 



